Decoherence in a Double-Dot Aharonov-Bohm Interferometer 



Bjorn Kubala, 1 David Roosen, 2 Michael Sindel, 3 Walter Hofstetter, 2 and Florian Marquardt 1 

1 Institute for Theoretical Physics, Universitat Erlangen-Niirnberg, Staudtstr. 7, 91058 Erlangen, Germany 
2 Institut fur Theoretische Physik, Goethe Universitat, 60438 Frankfurt /Mam, Germany 
3 Physics Department, ASC and CeNs, Ludwig- Maximilians- Universitat, 80333 Munich, Germany 

(Dated: November 15, 2010) 

Coherence in electronic interferometers is typically believed to be restored fully in the limit of small 
voltages, frequencies and temperatures. However, it is crucial to check this essentially perturbative 
argument by nonperturbative methods. Here, we use the numerical renormalization group to study 
ac transport and decoherence in an experimentally realizable model interferometer, a parallel double 
quantum dot coupled to a phonon mode. The model allows to clearly distinguish renormalization 
effects from decoherence. We discuss finite frequency transport and confirm the restoration of 
coherence in the dc limit. 
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Introduction.- Quantum coherence and its degrada- 
tion by interaction with the environment are crucial in- 
gredients in electronic transport through nanostructures. 
Generically, perturbation theory suggests that coherence 
is fully restored in the limit of small applied bias voltage, 
ac frequency, and temperature. In this limit, the parti- 
cles do not have enough energy to leave behind a "trace" 
in the environment. However, it is highly desirable to 
check whether such a statement remains valid beyond 
perturbation theory. In this work, we propose to employ 
the numerical renormalization group (NRG) [1] to study 
this question. That method can readily be applied to 
transport through interacting localized electronic levels, 
such as quantum dots and molecules. Investigating de- 
coherence in such systems offers the advantages of (some 
degree of) experimental control of the type and strength 
of decoherence mechanisms and, on the theoretical side, 
of a wealth of powerful, advanced methods. 

In discussing the effects of coupling a system to an 
environment true decoherence has to be carefully distin- 
guished from a mere renormalization of the system's char- 
acteristics due to the coupling. Electron-phonon (e-ph) 
coupling, for example, affects transport in several ways: 
Incoherent processes, in which phonons are absorbed or 
emitted, may increase or decrease the total transport 
through a single site [2], and coherent contributions to 
transport will also be renormalized. Considering an in- 
terferometer - the most obvious way to study decoherence 
in electronic transport - it thus turns out to be advanta- 
geous not to consider the overall change in the amplitude 
of the current or conductance as a measure of decoher- 
ence. Instead, we want to focus on the special case of 
destructive interference, i.e., the case of half a quantum 
of magnetic flux for an Aharonov-Bohm interferometer 
(ABI) (see Fig. [l]) . The advantage [3] is that in an inter- 
ferometer with perfectly symmetric arms renormalization 
processes will not lift the destructive interference, and 
any finite current is a true sign of decoherence. 

In this work, we consider a double-dot Aharonov-Bohm 



interferometer (see Fig. [T]), where electrons on the dots 
couple to a single phononic mode, representing the en- 
vironment. This scenario can be experimentally realized 
in tunneling transport through two quantum dots placed 
between two electron reservoirs [4] . On the one hand, one 
can then consider specifically engineered phononic envi- 
ronments, for instance, for quantum dots defined within 
free standing nanobeams [5] , or a quantum dot coupling 
to a cantilever [6]; on the other hand, the simple sce- 
nario can easily be generalized to model multiple phonon 
modes in standard heterostructures or molecular elec- 
tronics [7]. We will focus on the linear ac conductance, 
which is readily accessible by NRG and of great current 
interest in general [SJ d] . The nonperturbative NRG re- 
sults will be complemented by nonequilibrium Green's 
functions (NEGF) [JO] methods. 

The Model.- Our system is described by the Hamil- 
tonian H = H imp + J2 a =L/R( H ^s,a + H T , a ) with lead 

Hamiltonian H rcs . a = ^2k € kac\. a Cka, and spinless elec- 
trons on two degenerate dot levels, n — 1, 2 coupling to 
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Figure 1. A double-dot Aharonov-Bohm interferometer. De- 
coherence suppresses the amplitude of AB oscillations (dashed 
line). At destructive interference, (j>AB = $0/2, only inco- 
herent transport processes contribute in a symmetric setup. 
Mapping the Hamiltonian to (anti)symmetric combinations 
of the dot levels allows for a transparent exploitation of this 
cancellation of coherent transport. 
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a single phonon mode, 



[K a/3 (ui ac ) - K a p(0)} I (iw ac ) , where 



n=l,2 



where the level energy, e, is taken with respect to the 
Fermi level of the leads. Assuming an electron-phonon 
coupling of equal magnitude, but different signs for the 
two dots, the phonons effectively act as a kind of 'which- 
path' detector decohering the electrons. 

In a symmetric, destructive interferometer, i.e. for a 
tunneling Hamiltonian, H T , a = J2k,n=i,2 V a, n cl a d n + 
c.c, where V l /r, i = Vl,2 = -Vr,2 = V t it is deco- 
herence alone, which allows for transport, overcoming 
the destructive interference. This is seen by mapping 
the Hamiltonian to a new basis, a/s — {d\ =p obVv^, of 
(anti) symmetric combinations of the original dot levels. 
The impurity Hamiltonian then reads 



H imp = e (sU + a f a) + g(b + 6 f )(s f a + a f s) 



(2) 



with the (anti)symmetric level coupling to the left (right) 
lead only with tunneling matrix elements Vj,. s — V^ a = 
y/2Vt and Vr. s = Vr,, a = 0, see Fig. [I] It decouples into 
separated left and right parts for g — 0, whereas trans- 
port across the interferometer becomes possible due to 
a tunneling term between the two new levels, which is 
coupled to the phonon mode. While it is obvious that 
the noninteracting contribution vanishes, we will later 
show that destructive interference also cancels renormal- 
ized processes, reflected in the vanishing of certain classes 
of diagrams after performing the mapping. This particu- 
larly advantageous feature of only observing (in calcula- 
tion and experiment) incoherent transport processes of- 
fered by the setup considered here was not present in 
other theoretical studies of various interference effects in 
related setups [11] nor in studies of e-ph interaction for 
a single dot [I2HT5] . 

Techniques.- We will present two approaches to study 
this system. First, we employ non-equilibrium Green's 
function (NEGF) theory, where, based on the exact non- 
interacting result |16j , electron-phonon interaction is per- 
turbatively included using Keldysh diagrams. Particu- 
larly tailored to treat nonequilibrium situations, which 
can easily be realized in quantum dot systems by an ap- 
plied dc bias, the NEGF or Keldysh method [TU] has 
been widely used as a physically transparent tool to 
investigate e-ph interaction in such systems, often, as 
here, in a perturbative manner [12H14I or employing 
extensions of the polaron approach [T31 [TS], which al- 
lows tackling strong electron-phonon coupling for weakly 
tunnel-coupled systems. In this work, we focus on the 
linear ac- conductance G a p, giving the finite frequency 
current I a (uj &c ) = G a p{^ ac )Vp(uj ac ) flowing out of lead 
a as response to a small ac-excitation voltage Vp of 
frequency w ac applied to lead [T7]. It is connected 
to a (retarded) current-current correlator, <7a/3(w ac ) = 



^Q^(Wac) 



dte^i I a (t), I P (0) 



(3) 



To calculate this correlator, and, hence, the ac conduc- 
tance, we use a recently developed Meir-Wingrccn-typc 
formula for the linear ac conductance [TS], which al- 
lows straightforward diagrammatic calculations pertur- 
batively including electron-phonon interaction - com- 
plementing previous works building on a fully time- 
dependent Green's function formalism jTHj (see e.g., |20j . 
treating electron-electron interactions) and rate-equation 
approaches [2"T] . 

While the perturbative NEGF-approach is only valid 
for weak electron-phonon interaction, our second ap- 
proach, the numerical renormalization group method 
(NRG), invented by Wilson in the late 1970s [I], is espe- 
cially tailored for strongly interacting systems. Through- 
out the last decades numerous enhancements of the nu- 
merical renormalization group method have been intro- 
duced, ever improving the accuracy of correlation func- 
tions [2"2"H2"5] and opening the possibility to study quan- 
tum impurity systems coupled to a bosonic bath |26j . 
Here, we present to our best knowledge the first NRG 
calculation of ac conductance in a quantum dot system 
with an additional bosonic degree of freedom (cf. NRG 
studies on ac conductance with electron- electron interac- 
tion, for instance, in the Kondo regime |27j . and investi- 
gations of electron-phonon coupling effects on linear dc- 
conductance, spectral density and spin and charge sus- 
ceptibilities in [28, [29]). For a quantum impurity system 
with a fermionic bath and an additional bosonic mode 
coupled to the impurity directly, the numerical renor- 
malization group can be applied in a similar fashion as 
in the "pure" fermionic case without the bosonic mode. 
The mapping to the Wilson chain is not affected by the 
additional bosonic mode (25] [28] , since the bosons only 
enter in the very first step of the iterative diagonalization 
procedure. We implemented a NRG algorithm using the 
reduced density matrix [53] as well as the complete ba- 
sis of the underlying Fock space. Furthermore, we per- 
formed a z-trick averaging [30] over 32 slightly different 
discretizations of the conduction band to improve on the 
logarithmic resolution of energy (and thereby also fre- 
quency) around the Fermi energy inherent to NRG. In 
that manner, reliable numerical results up to frequencies 
of some ten percent of the bandwidth could be achieved 
with reasonably small numerical errors and artifacts. 

Results and interpretation.- To investigate the ef- 
fects of weak electron-phonon coupling on ac-transport, 
we employ NEGF, considering contributions to the ac- 
conductance from diagrams with up to a single phonon 
line, which gives contributions up to order 0{g 2 ). Using 
the mapping to (anti)symmetric dot levels in a diagram- 
matic representation for the ac conductance automati- 
cally exploits the main characteristic of transport in the 
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Figure 2. Mapping the Hamiltonian to (anti)symmetric com- 
binations of the dot levels: the mapped electron-phonon inter- 
action vertex (a) restricts the number and type of diagrams 
in a perturbative calculation of ac conductance (b), reflect- 
ing the destructive interference of all coherent contributions. 
(Solid lines denote electron Green's functions, wiggly lines 
stand for phonon propagators.) 



symmetric destructive ABI, namely the complete can- 
cellation of coherent transport processes due to destruc- 
tive interference. For instance, considering the ac con- 
ductance, ^Li?(w ac ), there are a zeroth-order and several 
second-order diagrams constructable from the original e- 
ph interaction vertex, Fig. [2] while no zeroth-order and 
only a single second-order diagram can be constructed 
from the mapped e-ph coupling vertex and the mapped 
lead-dot couplings. The nonexistence of those diagrams 
reflects the notion that in the destructive ABI the co- 
herent zero-order contribution, as well as contributions, 
where transport is "merely renormalized" , will interfere 
destructively and thus vanish. There are, however, dia- 
grammatic contributions for the "diagonal" conductances 
GLL/RRi^ac), describing electrons tunneling back and 
forth between the (anti)symmetric level and the (left) 
right lead, as the chemical potential of the lead is var- 
ied with frequency w ac . The noninteracting result thus 
found (see the lower right panels of Fig. 3]), which is 
also accessible by simpler techniques [§1 [T5]7reproduces 
for small frequencies the universal behavior of a meso- 
scopic RC-circuit with a quantized charge relaxation re- 
sistance [H E]. The second-order contribution to Gll 
constitutes a small correction to the noninteracting re- 
sult, which yields additional features at ui &c = e + ui^. 
Turning to the conductance across the interferometer, 
Glr, we find, as argued above, no zeroth order, and only 
a single, topologically distinct diagram in second order, 
namely an electronic loop with a crossing phonon line 
[see Fig. |2jb)]. One can now easily use Keldysh dia- 
grammatic rules to find the corresponding integral ex- 
pression for the ac conductance in terms of the elec- 
tronic Green's functions and the bare phononic propa- 
gators D^, A = (E' — w p h ± ir])^ 1 — (E' + w p h ± 277) 1 . 
Interestingly, there are contributions from the <5-peak at 
the phonon energy oj p h as well as from the principal- value 
part. For the dissipative real part of the conductance, 
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Figure 3. ac conductance through the destructive ABI. In the 
limit of weak electron-phonon coupling NRG and perturbative 
calculations match. For stronger e-ph coupling, g/uj p h ~ 1, 
additional sidepeaks at oj ac = e + nuip^ (n = 1, 2, . . .) ap- 
pear, and the first resonance is shifted to lower frequency. 
Glr(uh.c = 0) = is found also for strong e-ph coupling. Pa- 
rameters are w ph = 0.05, e = 0.025, T = 2ix(p L + pn)\ V t \ 2 = 
0.02 in units of the bandwidth, conductances in units of e 2 /h, 
and we consider the zero temperature limit and no dc-bias. 



Re<7i_R, (see left panel of Fig. |3j), we find that contribu- 
tions from the latter yield a resonant peak at w ac = e. 
Around ui ac = e + w p h, where the excitation frequency 
is in resonance with the dot level position after the ex- 
citation of a phonon, there is a smaller negative peak, 
stemming from the ^-contributions. Corresponding fea- 
tures at the same frequencies are found for the imaginary 
part, lmQ LR . 

Investigating strong electron-phonon coupling, compar- 
isons of NRG calculations to the perturbative results 
are shown in Fig. [3| For weak e-ph coupling, there 
is good qualitative and also quantitative agreement be- 
tween NRG and the perturbative results: In the limit 
g \ 0, we reach the noninteracting results for Gll (thick 
black lines in the two lower right panels of Fig. |3| , while 
for Glr we recover the perturbative, 0(g 2 ), result derived 
above. The left panel of Fig. [3] also demonstrates the 
main effects of stronger e-ph coupling on Re Glr- Most 
importantly, one finds additional resonances at sidebands 
of the level energy, w ac = e + riw p h (n = 1, 2, . . .), which 
stem from processes involving several phonons. Fur- 
thermore, we observe a shift of the first resonance from 
^ac ~ e to smaller frequencies reminiscent of a polaron 
shift. 

In the following, we want to highlight two particular 
remarkable features of our results at low frequency. The 
first, observable already in the weak coupling limit, is the 
fact, that ac transport is already possible for excitation 
frequencies smaller than the phonon energy, w ac < lj p ^ . 
Let us first emphasize, that this is in striking contrast 
to the case of finite-bias dc transport, where in the 
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weak coupling limit there is a sharp threshold for in- 
coherent transport processes becoming possible only for 
eVdc > Wph- Using the mapping to (anti)symmetric dot 
states the dc current through a destructive ABI can be 
shown to contain only incoherent contributions. These 
contributions are characterized by explicit appearance 
(cf. Eq. 66 in Ref. 12 ) of self-energies, £§:>P h , due to 
emission or absorption of phonons. In lowest order in the 
e-ph coupling (and for zero temperature) transport takes 
place through two "transport channels" at uj = e, e+Wph 
within a transport window, fi^ > u> > [1r + u; p h . In par- 
ticular, the resulting I-V curve has a sharp gap at low 
bias, eVdc < ^ph, as any incoherent transport process 
necessarily involves the emission of a phonon, whose en- 
ergy has to be provided by the applied bias. The naive 
expectation, that in ac transport the frequency of the 
infinitesimal ac-excitation voltage providing energy for 
emission of phonons acts in a similar manner does not 
hold. In fact, subthreshold ac transport apparently is 
possible without the real emission of phonons. Con- 
sidering the Hamiltonian of the destructive ABI in the 
mapped base, we see that virtual tunneling back and 
forth processes between symmetric and antisymmetric 
levels with emission and absorption of phonons will yield 
an effective charge-charge interaction between the two 
sites. Such interaction then leads to finite ac conduc- 
tance, but vanishing dc conductance, just like the elec- 
trostatic interaction in a plate capacitance. 

While both of these arguments explain the vanishing 
of conductance in the dc limit in lowest order in the e- 
ph coupling, note as a second, most important feature 
of the results in Fig. [3j that our NRG results indicate 
that this holds for strong coupling also. Indeed, while 
in a lowest-order process phonons can only be excited 
exactly at their bare frequency w p h , considering higher 
orders, one may naturally ask about the effects of the 
broadening of the phonon propagator due to coupling 
to the electrons. Some (though not all) higher-order ef- 
fects can thus be captured by calculating the electronic 
self-energies with a phonon propagator dressed with elec- 
tronic polarization bubbles. In the dc-case, one finds that 
electrons can indeed absorb an energy up to eV from the 
phononic bath broadened by coupling to the electrons, in- 
stead of needing to emit w p h as in lowest order. Nonethe- 
less, the weight of such contributions is small leading to 
a current dependence I cx (eV) 3 for small bias voltage 
for the considered higher-order contributions in the e- 
ph coupling and consequently a vanishing linear conduc- 
tance at zero frequency. This agrees with the dc limit 
of our NRG results for the ac conductance, which indi- 
cates as discussed that all higher-order contributions to 
the linear zero-temperature conductance vanish. We thus 
confirm the naive perturbative argument, that there are 
no incoherent transport processes at vanishing frequency, 
temperature and transport voltage. 



Conclusions.- We have studied decoherence and the 
effects of electron-phonon coupling on ac transport in 
a double-dot interferometer. A diagrammatic calcula- 
tion reveals that renormalization processes can be clearly 
distinguished from decoherence for the case of destruc- 
tive interference, and it explains low frequency transport 
as a consequence of phonon-mediated charge-charge in- 
teractions in a mapped system. NRG calculations for 
stronger coupling show characteristic features due to 
multi-phonon processes. With increasing coupling, the 
phonon mode broadens and effectively yields a contin- 
uous environment spectrum. Nevertheless, the nonper- 
turbative NRG result confirms that coherence is fully re- 
stored in the zero-frequency limit of the linear conduc- 
tance at zero temperature. 
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